Richtmyer-Meshkov instability and solid 4 He melting driven by acoustic pulse 
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Recent experiments have shown remarkable dynamics of solid 4 He melting and growth, driven 
by the normal incidence of an acoustic pulse on the solid-liquid interface. The theory of solid 
growth/melting, driven by the radiation pressure of the acoustic pulse, accounts well for the temper- 
ature dependence of the measured data. There is however an observed source of extra, temperature- 
independent, melting. We here propose that this extra melting is due to solid-liquid mixing (and 
consequent melting) at the interface, in a process similar to the Richtmyer-Meshkov instability: Ini- 
tial undulations of the rough interface, grow when accelerated by the acoustic pressure oscillations. 
This model predicts a temperature-independent extra melting and its dependence on the acoustic 
power, which is in agreement with the measured data. 

PACS numbers: 67.80.-s,68.08.-p,43.35.+d,47.20.k 



Recent experiments [l| have shown that the radiation 
pressure P r of an acoustic pulse, incident on the solid- 
superfluid interface of 4 He, can induce the local melting 
or growth of the solid. The general temperature depen- 
dence of the pulse-induced melting/growth is accounted 
for using the concept of acoustic radiation pressure. Nev- 
ertheless, a large deviation was found, which seems to be 
temperature independent. It is this extra melting which 
we shall describe in this paper. 

We begin with the analysis of the experimental data 
in terms of a linear growth/melting coefficient (Eq.1-6 of 
0). The amplitude of melting/growth is linearly related 
to the radiation pressure P r through the (temperature 
dependent) growth coefficient K(T) [ij 
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where r = 1msec is the overall duration of the acoustic 
pulse. The radiation pressure is itself temperature de- 



pendent: P r = E 
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where the en- 



ergy density of the incident wave is E — I/c\, the power 
density / = 144W/m 2 , and the reflection coefficient is 



Z 2 - Zi 



R = 



z\z 2 p\ 



/'■2 



z± + z 2 + 



z 1 z 2 p 1 K(T) 



Pl~P2 
P1P2 



(2) 



The growth coefficient K(T) is 



where zi, 2 — p\, 2 c\ y2 . 
given by 

K (T)- 1 = AT 4 + Be'^' (3) 

(with A = 3 x 10- 2 ,B = 2.4 x 10 3 ,A = 7.2K 2]). In these 
expressions pi. 2 , ci : 2 are the density and sound velocities 
of the two phases on either side of the interface, with the 
index 1 corresponding to the phase on the side of the inci- 
dent wave. In the experiment [IJ we have pi = 0.17gr/cc, 
ci = 365m/sec and p s = 0.19gr/cc, c s — 475m/sec. In 
Fig.l we compare the experimental data with the cal- 
culated growth/melting amplitude h, given by Eq.JIJ. 



The overall agreement is good (dotted lines in Fig.l), 
but there is clearly an additional source of melting which 
appears for pulses incident from both the solid and liq- 
uid sides. This extra melting is not caused by the sound 
heating of the solid |3J. 

The contribution from the acoustic pressure leads to 
growth/melting that proceeds at a constant rate for the 
duration of the applied pulse Q . In addition there is the 
inertial acceleration of the solid- liquid interface by the in- 
cident pressure oscillations. This motion should average 
to zero if no instability or mixing occurs at the interface. 
It is only due to mixing at the interface that this mecha- 
nism can contribute on average to the position of the in- 
terface at the end of the pulse. We assume here that the 
shape changes of the interface due to Richtmyer-Meshkov 
(RM)-mixing Q do not change the average acoustic pres- 
sure and therefore can be approximately decoupled from 
each other. We therefore add the contribution from the 
acoustic pressure to the mixing contribution due to the 
RM-mixing. A complete description has to take into ac- 
count possible interference effects of one mechanism with 
the other. 

Let us now summarize the main results of this paper. 
We propose that the measured extra melting is the result 
of dynamic solid-liquid mixing at the interface. This mix- 
ing is due to a RM- instability Q|, that typically occurs 
when two liquids/gases of different densities are accel- 
erated one against the other by an incident shock-wave, 
causing initial perturbations of the interface to grow. In 
our case the interface is between two phases of the same 
material, but of different density, which is repeatedly ac- 
celerated to and fro by a series of rapid pressure oscilla- 
tions. Under these conditions, and a rough initial inter- 
face, solid-liquid mixing should occur due to a RM-type 
instability |4|, since we note that in the limit of weak 
shock- waves these behave as acoustic pulses [5} . Using 
the linear growth theory of the RM instability we predict 
a temperature-independent melting, with the correct de- 
pendence on the incident acoustic power. With a single 
free parameter we get good quantitative account of the 
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measured extra melting. 

A RM instability will induce a layer of mixing at the in- 
terface between two materials of different densities (here 
the solid-liquid interface). The mixing occurs when ini- 
tial perturbations of the interface between two liquids (or 
gases) grow after the interface is accelerated by a shock- 
wave Q (Fig. 2). The linear single- mode theory describes 
the inertial growth of an initial interface undulation, due 
to an impacting shock-wave, gives 

— = kaovA =>- clrm = »o + kaovAr (4) 

where k is the wavevector of the interface undulation of 
initial amplitude ao, v is the interface velocity due to the 
acceleration and A = (p 2 — P\)/{p2 + Pi) is the Atwood 
number. The sign of A relates to the phase-reversal pro- 
cess (Fig. 2b) and is irrelevant here. Note that Eq.(0J is 
temperature-independent. Since in our case the solid- 
liquid interface is continuously accelerated to and fro by 
the pressure oscillations |lOj, we therefore interpret the 
equation for the rate of growth of surface instabilities 
as giving the rate at which a continuous flux of solid is 
locally mixed with the liquid and melts. 

Furthermore, unlike the usual RM scenario, we have 
a solid crystal on one side of the interface. Due to 
the rapid melting-solidification at the rough (i.e. un- 
faceted) solid-superfluid interface, this interface deforms 
similar to that between two liquids pj. On the other 
hand, dissipation and surface tension tend to damp out 
any inertial growth of instabilities. The solid-liquid sur- 
face tension, is a ~ 0.2erg/cm 2 for 4 Hc 8], and acts 
to stabilize and damp the inertial growth of initial per- 
turbations of the RM-type [jj. The regime where the 
surface tension dominates occurs whenever the damping 
rate u> = k 3 ^ 2 yja/ {fix + p 2 ) is larger than the growth 
rate given in Eq.Q, i.e. for k > (a vA) 2 (pi + p 2 )/a. 
Using the relevant values (estimating ao < 10 A), we find 
that in the experiments the surface tension dominates. 
This means that following a single pressure oscillation, 
the amplitude growth described in Q would have been 
quickly damped out, and result in negligible solid-liquid 
mixing. 

Nevertheless, in the experiment p], there is not a sin- 
gle pressure oscillation followed by inertial growth, but 
rather a long pulse of continuous acoustic waves. Each 
pressure oscillation (of ~ 0.1/isec duration, frequency 
9.9MHz 0) will cause the solid- liquid interface to ac- 
celerate, and existing undulations to grow in amplitude. 
The solid-liquid interface is therefore continuously ac- 
celerated back and forth, throughout the pulse dura- 
tion. The damping due to the surface tension is therefore 
not relevant, and we proceed to estimate the amount of 
solid-liquid mixing due to this continuously driven RM- 
instability, using the simplest linear growth theory Q 
(Eq. 2J), over the entire pulse duration r. In this case we 
treat the growth rate equation (0} as a measure of the 
rate at which solid material is continuously mixed into 
the liquid, in a manner similar to the (highly exagger- 



ated) "spike" in Fig. 2. This material then melts as it 
finds itself in a region which is liquid in equilibrium. The 
accumulated amplitude cihm represents therefore the to- 
tal thickness of solid that was mixed (and melted) by 
this process over the duration of the acoustic pulse. In 
particular, it does not correspond to the final amplitude 
of individual interface undulations. 

The interface velocity v, due to the oscillating pressure 
of amplitude SP = J ' Ep\^ 2 c\ 2 (depending on the direc- 
tion of the incident wave), is given in terms of the local 
material velocities at the interface, after the passing of 
the pressure wave 0: v = (piv m ,i - P2V m a)/{pi ~ Pi)- 
The local material velocities following a pressure oscilla- 
tion are given by ^(| v m = dP/ z i,2- For power density 
I = 144W/m 2 we therefore get: v ~ 0.023m/sec. 

In Fig.l we compare the calculated combined melt- 
ing/growth h — clrm HI Hi with the experimental data. 
We find that a good agreement 01 is achieved when 
we use an interface roughness aspect-ratio of kao ~ 20 
in Eq.Q. For a rough solid-liquid interface, which is 
the case in the experiments done so far 0, Q , we would 
naively expect kao ~ 1. This discrepancy indicates that 
the linear growth theory is not entirely adequate, espe- 
cially as it usually deals with a single pressure pulse, 
while in our case the interface is continuously driven. 
The use of a constant amplitude ao in the rate equation 
Q is a rough approximation, as the amplitude of surface 
undulations grows between each pressure oscillation. We 
therefore note kao as the free parameter in our model, 
and takes into account our ignorance of the details of the 
solid-superfluid interface dynamics. Note that the agree- 
ment is better for the case of melting with an acoustic 
pulse from the liquid side (Fig.l). 

Our model further predicts that olrm °c Sv m cx vJ, 
while the radiation-pressure driven growth predicts a 
linear dependence on the acoustic power. In the regime 
where the melting is dominated by the RM-instability, 
i.e. at relatively high temperatures, we therefore expect 
to find a RM cx / behavior. The experimental data at 
1.2K [3( does indeed show such a power law (Fig. 3). 

To conclude, we have analyzed recent experiments 
where the solid-superfluid interface is accelerated by an 
acoustic-pulse, and a large source of unexplained melting 
was observed. It was suggested that some mechanism of 
energy dissipation at the interface is responsible for this 
unexplained melting 1]. We proposed here a possible 
mechanism for such a dissipation; local solid-liquid mix- 
ing (and melting) occurs due to a dynamic instability 
of the interface, of the RM-type. The proposed model 
is in good agreement with the experimental data, and 
describes a specific mechanism by which the accelerated 
interface melts the solid. 

It will be interesting to study this instability using sim- 
ulations taking into account both the material flows 
and the solid-liquid melting/solidification dynamics. Ex- 
periments on faceted crystals, with an atomically smooth 
interface, should show less melting due to the RM-like 



3 



process described here. It could also be useful to repeat 
these experiments with a controlled initial undulation on 
the solid-liquid interface, of macroscopic size, and visu- 
alize its evolution when accelerated by an acoustic pulse. 
This initial undulation can be formed by a standing cap- 
illary wave at the interface fill ]. 
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FIG. 1: The measured melting/growth due to an acoustic 
pulse from the solid/liquid side £jj (circles & squares re- 
spectively). We plot the calculated radiation- pressure driven 
growth/melt h due to the acoustic pressure (Eq0, using 
A = 7.2K (dotted lines) and A = 7.5K (solid lines) in (Eq^J. 
The dashed lines show the result of subtracting the RM- 
amplitude orm (since it represents extra melting) (Eq^J from 
the calculated linear growth/melt ( Eq|S| using A = 7.5K, 
solid lines), using kao — 20. 
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FIG. 2: An illustration of the evolution of a Richtmyer- 
Meshkov instability |(J and the consequent material mixing 
(bottom line). An initial sinusoidal undulation on the in- 
terface (top line) is accelerated by an impacting shock-wave 
(horizontal line with the attached arrow), from either the (a) 
light (low density) or (b) heavy (high density) side. The in- 
ertial mixing shown here ("spike" in the bottom line), after 
an acceleration by a single pressure pulse, corresponds to the 
case of negligible surface tension, unlike the solid-superfluid 
4 He case. 
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FIG. 3: Measured melting depth as a function of the acoustic 
power (both normalized by their largest measured values) at 
T=1.2K (squares 3D- The straight line is a fit to illustrate 
the a 2 RM oc I behavior, as predicted by our model. 



